J/tfj azimuthal anisotropy relative to the reaction 
plane in Pb-Pb collisions at 158 GeV per nucleon 
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Abstract 

The J/t/> azimuthal distribution relative to the reaction plane has been measured 
by the NA50 experiment in Pb-Pb collisions at 158 GeV/nucleon. Various physical 
mechanisms related to charmonium dissociation in the medium created in the heavy 
ion collision are expected to introduce an anisotropy in the azimuthal distribution of 
the observed i/tp mesons at SPS energies. Hence, the measurement of J/i/j elliptic 
anisotropy, quantified by the Fourier coefficient V2 of the J/ip azimuthal distribution 
relative to the reaction plane, is an important tool to constrain theoretical models 
aimed at explaining the anomalous J/1/1 suppression observed in Pb-Pb collisions. 
We present the measured J/?/> yields in different bins of azimuthal angle relative to 
the reaction plane, as well as the resulting values of the Fourier coefficient V2 as 
a function of the collision centrality and of the J/ip transverse momentum. The 
reaction plane has been estimated from the azimuthal distribution of the neutral 
transverse energy detected in an electromagnetic calorimeter. The analysis has been 
performed on a data sample of about 100 000 events, distributed in five centrality 
or pt sub-samples. The extracted V2 values are significantly larger than zero for 
non-central collisions and are seen to increase with p-r. 



1 Introduction 

Charmonium production and suppression is one of the most powerful probes for 
a phase transition to deconfined matter in heavy-ion collisions at the energies of 
the CERN SPS. In particular, the 3/tp suppression in proton-nucleus and nucleus- 
nucleus reactions has been intensively studied in the last 2 decades by the NA38 
and NA50 experiments. The J/1/1 suppression observed in proton-nucleus reactions 
is understood as due to absorption of charmonium states on ordinary nuclear matter 
with <j a bs = 4.2 ± 0.5 mb Q~]. The J/i/>/Drell-Yan ratio measured in S-U and 
peripheral Pb-Pb collisions results to be in agreement with the expectation from 
ordinary nuclear absorption as measured in p-A reactions, while an anomalous extra 
suppression is present in semi-central and central Pb-Pb collisions [2J. 



1 



Additional insight into charmonium suppression mechanisms can be obtained 
from the anisotropy of the overlap region of the projectile and target nuclei in col- 
lisions with impact parameter b> |3l |4| |]. The initial geometrical anisotropy 
gives rise to an observable anisotropy in particle distributions if the created system 
is interacting strongly enough to thermalize at an early stage and develop collec- 
tive motion (flow). Hence, anisotropic transverse flow should be observed for J ftp's 
formed by c-c recombination if the charm quarks have undergone strong enough re- 
scatterings leading them to thermalize in the first stages of the system evolution. 
J/ip flow is however not expected to be established at SPS energies where early 
charm thermalization is unlikely and cc recombination is negligible. Nevertheless, 
other mechanisms related to cc absorption in the medium created in the collision, 
essentially cc dissociation by the hard gluons present in the deconfined phase [4] [5] 
and by co- moving hadrons [3] are predicted as possible sources of 3 ftp anisotropy 
already at SPS energies. 

The azimuthal anisotropy is usually quantified from the coefficients of the Fourier 
series describing the particle azimuthal distribution: 



where ^rp is the reaction plane angle defined by the impact parameter vector in the 
transverse plane. An azimuthal dependent J ftp absorption pattern determined by the 
anisotropic geometrical shape of the nuclear overlap region is expected to give rise to 
a measurable second harmonic coefficient V2, which describes an elliptic anisotropy. 
It is anticipated [4j [5] that elliptic anisotropy due to J/ip dissociation by gluons 
resulting from the formation of a deconfined medium should vanish for peripheral 
collisions (where the critical temperature is not attained) and for central collisions 
(because of the isotropic geometry of the overlap region), showing a sudden onset 
in correspondence of the phase transition and a maximum for semi-central collisions 
(£ T ~70-80 GeV in g] or iV part -200-220 in 0). 

2 Experimental setup, data selection and reac- 
tion plane estimation 

The NA50 apparatus consists of a muon spectrometer equipped with three detectors 
to measure centrality-related observables on an event-by-event basis and specific 
devices for beam tagging and interaction vertex identification. A detailed description 
of the detectors can be found in [6j. Minimal details are given here for specific 
detectors relevant for the present analysis. 

Anisotropy studies have been done on the sample of about 100,000 J/tp's col- 
lected by experiment NA50 in year 2000 with the SPS Pb beam at 158 GeV/nucleon 
(see [2]). The J ftp is detected via its fi + \x~ decay in the pseudo-rapidity range 2.7 < 
?7iab< 3.9. The analysis is performed in the dimuon kinematic domain < j/ cm < 1 
and —0.5 < cos(#cs) < 0.5, where y cm is the rapidity in the center-of-mass system 
and 8cs is the polar decay angle of the muons in the Collins-Soper reference frame. 

The study of the centrality dependence of the J ftp anisotropy is based on 5 bins of 
the neutral transverse energy (Et) as measured by an electromagnetic calorimeter 
(see fig. [1]). This calorimeter is made up of lead and scintillating fibers and it 
measures event-by-event the transverse energy carried by neutral particles produced 
in the interaction (mostly due to it — > 77 and to direct 7) in the pseudo-rapidity 
window 1.1 < 77i a b< 2.3. The Et limits for the 5 bins are reported in table[T]together 
with the average values of impact parameter and number of participants estimated 
by means of a Glauber calculation including the resolution of the calorimeter. 

The reaction plane is estimated making use of the azimuthal segmentation in 
six azimuthal sectors (sextants) of the electromagnetic calorimeter, as represented 
in fig [T] where it can also be seen that each sextant is further subdivided into four 
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Tabic 1: Et limits, average values of impact parameter and number of participants for 
the ccntrality bins used in this analysis. 




Figure 1: Geometrical segmentation of the electromagnetic calorimeter. 



radial rings, each of them covering a different pseudo-rapidity range. The event 
plane angle ^ n (estimator of the unknown reaction plane angle ^ rp) is given by: 



— —tan 

n 



Eti #t sin(n$Q 
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where n is the considered Fourier harmonic, the neutral transverse energy mea- 
sured in sextant i, and $i the azimuthal angle defined by the center of sextant i 
(see fig. [TJ . The weighting coefficients w l are introduced to make the event plane 
distribution isotropic and are defined as < E^ ot > /(6 < E l T >). Their values range 
between 0.994 and 1.012, resulting in a very small event-plane flattening correction. 
The event plane ^ has been used to calculate the elliptic anisotropyQ It is com- 
puted from the ir° azimuthal distribution in the backward rapidity region, where, at 
SPS energies, pions show positive «2 [H [H E2 an d therefore it is directed in-plane 
(i.e. parallel to the reaction plane). 

The event plane resolution (expressed as < cos [2 (^2 — ^rp)] > ) has been 
estimated in two independent ways. The results are shown in fig. [2b. The first 
technique is based on Monte Carlo simulations of the detector response taking as 
input the value of V2 measured by the calorimeter [10]. The sextant to sextant 
fluctuations are tuned to reproduce the experimentally observed distribution of the 
quantity: 

03 = : — 3^ (a) 

7rsinc^f 

6 

with sine x = (sinx-)/x. For symmetry reasons, 63 is sensitive only to statistical 
fluctuations and not to azimuthal anisotropies [TD]. The gray band represents the 



1 The event plane could also be used, but the resolution is worse [7|. 
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systematic error coming from the systematic uncertainty on vi due to the presence 
of non-flow correlations. The second technique makes use of the ring segmentation 
of the calorimeter to define two sub-events. The resolution is extracted from the 
angular correlation between the event plane angles of the two sub-events. A scheme 
of the sub-event definition is shown in fig. [2^,: ring 2 has been removed from the 
analysis in order to have a rapidity gap limiting non-flow correlations between the 
two sub-events, while the alternate pattern of rings and sextants is dictated by the 
need of having the two sub-events equally populated. The energy collected by the 
excluded ring 2 is then accounted for when using the formulas from [7] to extrapolate 
the measured resolution of the sub-event plane to the resolution of the event plane of 
the full calorimeter. As it can be seen in fig. Hp the resolutions extracted with the two 
methods agree within the systematic uncertainties, so the Monte Carlo estimation 
together with its systematic error bar is used to calculate the J/t/j elliptic anisotropy. 
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(a) (b) 

Figure 2: (a) Sketch of the geometrical configuration used for sub-event definition, (b) 
Second harmonic event plane resolution with from Monte Carlo simulations and sub-events 
technique. 



3 Analysis and results 

The J/ip azimuthal distribution relative to the reaction plane is extracted using two 
different analysis schemes. The first one extracts the number of J/i/i's in different 
intervals of azimuthal angle relative to the event plane. The second one computes 
the Fourier coefficient vi from the J/ip azimuthal distribution. For each analysis, 
specific methods are implemented in order to separate the J/ip signal from the other 
dimuon sources under the J/ip mass pealfl 

3.1 Number of J/^'s in bins of azimuthal angle 

Two different analysis methods have been developed to extract the number of J ftp's 
in two wide bins of azimuthal angle relative to the event plane (A$2 = $jim« — ^2 
where $d, mu is the dimuon azimuthal angle and *2 the second harmonic event plane 
from eq.[5J). 

The first method consists in fitting the mass spectra of opposite-sign dimuon 
sub-samples in bins of centrality (Et) and dimuon azimuthal angle relative to the 
measured event plane (A$2)- The mass spectra above 2.5 GeV/c 2 (see fig. O-left) are 
fitted to the four signal contributions (namely 3/ijj, ip', Drell-Yan and open charm) 
with shapes determined from detailed Monte Carlo simulations of the NA50 appara- 
tus. The combinatorial background is evaluated from the like sign pairs (see [2] for 

2 About 90% of the reconstructed opposite-sign dimuons in the mass range 2.9 < M < 3.3 GeV/c 2 and 
in the kinematic domain defined above originate from J/ip decays. 
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Figure 3: Left: Example of fit to the mass spectrum. Right: Et spectra (not 

corrected for ccntrality-depcndent inefficiencies) of \x~ in 2.9 < M < 3.3 GcV/c 2 . 



details) . This analysis method is limited by the low statistics of high- mass Drell-Yan 
dimuons which are crucial to fix the Drell-Yan contribution in the fitting procedure. 
Hence, it is not possible to divide the full sample of dimuon events in more than 10 
bins (5 centrality x 2 A<I>2 intervals) 

The second method consists in building, for each A$2 bin, the Et spectrum 
of all the [i + in the mass range 2.9 < M < 3.3 GeV/c 2 and then subtracting 
the spectra of the different sources of background. The Et spectra of the various 
dimuon contributions in 2.9 < M < 3.3 GeV/c 2 are shown in fig. [3} right. The com- 
binatorial background is extracted from like-sign muon pairs. The DY spectrum 
is estimated from fj, + fi~ in the mass range 4.2<M<7.0 GeV/c 2 and rescaled, via 
Monte Carlo simulations, to 2.9 < M < 3.3 GeV/c 2 . The DD yield is estimated 
from opposite sign dimuons in 2.1<M<2.7 GeV/c 2 after combinatorial background 
and DY subtraction and rescaled to the J/if> mass range. The dimuons from i/j' de- 
cay in 2.9 < M < 3.3 GeV/c 2 are negligible. The underlying assumption is that the 
Et spectra of the involved physical processes do not depend on the invariant mass 
range used for their determination within the range under study. This "counting" 
method allows for a larger number of Et bins, thus providing a better insight on 
the centrality dependence of the possible anisotropy. It should be noted that no cor- 
rection for centrality-dependent inefficiencies is applied to the Et spectra because 
in the following analyses only relative numbers of J/^'s in different A<I>2 bins are 
considered. 

The presence of an azimuthal dependent J/V> absorption is expected to result in 
a different number of particles emitted parallel (in-plane) and orthogonal (out-of- 
plane) to the reaction plane as a consequence of the geometrical shape of the overlap 
region of the colliding nuclei. The elliptic anisotropy is therefore quantified starting 
from the numbers Nin and Nout of J/tp's observed in two cones with an opening 
angle of 90° centered respectively at A$2=0° (in-plane) and at 90° (out-of-plane), 
see fig. |4] So: 

NlN = / ,, A . r d A$ 2 + / d A$ 2 

The anisotropy is then quantified as the ratio (Nin — Nout) / (Nin + Nout)- 
A positive anisotropy comes from a larger number of J/if>'s observed in plane than 
out-of-plane. If only a second (elliptic) harmonic is present, i.e. dN / d&dimu oc 
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Figure 4: Anisotropy vs. -Et as extracted from the number of J/tjj's in A(p„ bins centered 
in-plane and out-of-plane (right). Error bars represent statistical errors on the measured 
points. 



1 + 2v 2 cos[2($ dimll - *flp)], then: 

N IN - Nout 4 . . 

m — Tat = Z V2 W 

JV/AT + IS OUT TT 

It should be noted that the resolution of the event plane has not been taken into 
account in the calculation of the anisotropy from Nin and Nout and therefore the 
comparison with the values of the Fourier coefficient v 2 reported in the next section 
is not straightforward. The results for the elliptic anisotropy as a function of Et are 
shown in fig. [4] The two analyses agree in indicating on average a small excess of 
3 ftp's emitted in-plane (positive anisotropy). The largest signal is observed in the 
centrality bin 70<E T < 90 GeV. 



3.2 Fourier coefficient v 2 of J/^'s 

The second coefficient of the Fourier expansion is given by: v 2 =< cos[2($ £ j im „ — ^ rp)] > 
where the average is performed over events in a given centrality (or p-r) bin. Since the 
reaction plane (^rp) is unknown, the event plane (calculated from neutral transverse 
energy anisotropy) has to be used instead, obtaining v' 2 =< cos[2($* mtl — ^2)] >• 
The quantity v 2 should then be corrected for the event plane resolution [7], obtaining 
V2 = v' 2 / < cos[2(* 2 - Vrp)] >■ 

Two different analysis methods have been used to subtract the background and 
extract the values of 3/tp elliptic anisotropy v%. The first estimation is obtained 
from the average of the cos[2($<jj mn — $2)] distributions of fi + fi~ in the mass range 
2.9 < M < 3.3 GeV/c 2 after subtracting the background contributions with the same 
"counting" procedure described above. The cos[2($di mtl — 9 2)] distribution of com- 
binatorial background is estimated from like-sign muon pairs, while the ones of DY 
and DD are extracted from different fi~ mass intervals and rescaled to the J/if> 
mass range under the assumption that they do not depend on the invariant mass 
range considered. 

A second evaluation of v 2 in the 5 Et bins has been obtained from the number 
of J ftp's extracted with the "counting" method in 8 bins of azimuthal angle relative 
to the event plane. The coefficient v 2 is obtained by fitting the resulting number of 
J ftp's in bins of A"I>2 with the function 

N J/ ^(A$ 2 ) =K[l + 2-v' 2 cos(2A$ 2 )] (7) 

where A$2 = &dimu — ^2 and the free parameters of the fit are K and v' 2 . Afterward, 
v 2 is obtained by applying to v' 2 the correction factor for the event plane resolution. 
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Figure 5: Fourier coefficient i>2 vs. Et (left) and 3/ip px(right). Error bars represent 
statistical errors, the gray band is the systematic error coming from the estimation of the 
event plane resolution. 



The results for V2 vs. Et are shown in fig. [5]- left. The two analysis methods 
are in remarkable agreement and show positive values of v%, confirming the excess 
of J/^'s exiting in-plane. A maximum V2 is observed for the bin 70<i?T<90 GeV, 
corresponding to N par t ~ 270 and (b}=4.8 fm. The error bars represent the statis- 
tical error on the measurement of the 3/ip anisotropy, while the gray bands are the 
systematic errors coming from the uncertainty on the estimation of the event plane 
resolution. The analysis based on the cos[2($dimu — ^2)] spectra has been applied 
also in bins of J/-0 transverse momentum. The obtained results for V2 as a function 
of pt (centrality integrated) are shown in [fright. The J/ip V2 shows an increasing 
trend with increasing pr. 



4 Conclusions 

J/ip elliptic anisotropy relative to the reaction plane has been measured by NA50 
from a data sample of 100000 J/i/i's produced in Pb-Pb collisions at 158 GeV/nucleon 
(y/s = 17.2 GeV). The anisotropy has been quantified both from the normalized 
difference between the number of J ftp's emitted in plane and out-of-plane and from 
the Fourier coefficient (V2) which describes an elliptic anisotropy. These quantities 
have been measured as a function of collision centrality (defined by the neutral 
transverse energy Et produced in the collision) and as a function of 3/ip transverse 
momentum. A positive V2 is measured: more J/^'s are observed in-plane than out- 
of-plane. The largest anisotropy is observed in the centrality bin with N par t ~ 270 
and (b)=4.8 fm. The elliptic anisotropy is observed to increase with increasing J/ip 
Pt- 
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